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Abstract 



We introduce a self-avoiding walk model for which end-effects are completely eliminated. We enumerate the 
number of these walks for various lattices in dimensions two and three, and use these enumerations to study the 
properties of this model. We find that endless self-avoiding walks have the same connective constant as self- 
avoiding walks, and the same Flory exponent v. However, there is no power law correction to the exponential 
number growth for this new model, i.e. the critical exponent 7=1 exactly in any dimension. In addition, the 
number growth has no analytic corrections to scaling, and we have convincing numerical evidence to support the 
conjecture that the amplitude for the number growth is a universal quantity. The technique by which end-effects 
are eliminated may be generalized to other models of polymers such as interacting self-avoiding walks. 

Keywords self-avoiding walk; critical exponent; enumeration; critical amplitude 



Self-avoiding walks (SAW) on a regular lattice are an important model in statistical mechanics with a long his- 
tory [1], An n-step SAW is a map cj from the integers {0, 1, • • • , n} to sites on the lattice, with cj(0) conventionally 
at the origin, and \u>(i+l) — u>(i)\ = 1, and co(i) 7^ U!(j)Vi 7^ j. The SAW model is the N — > limit of the TV-vector 
model, which includes the Ising model as the N = 1 case, and as such it is an important model for the theoretical 
understanding of critical phenomena [2]. Self-avoiding walks are in the same universality class as polymers in a 
good solvent, and thus the model also has an important physical interpretation. Much is known about self-avoiding 
walks, particularly in dimensions d > 4 [3, 4], and considerable progress has been made in understanding self- 
avoiding walks in two dimensions [5, 6, 7], but an exact solution remains elusive for any dimension bar the trivial 
case d = 1. 

The most important quantities which characterize SAW are the number of SAW of length n, c„, and measures 
of the size of the walk, such as the square end-to-end distance, r n . In the following definitions we assume that 
summations are performed over the set of self-avoiding walks of length n: 



We abuse notation by allowing to measure the number of steps of a walk oj, while |w(n)| 2 is the square of the 
Euclidean norm of the site co(n). 

The asymptotic behavior of c„ and r„ is characterized by certain quantities of fundamental interest, namely the 
connective constant, [i, which is lattice dependent, and the critical exponents 7 and v which are universal. We have 
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where we define the mean square end-to-end distance as p n = r n /c n . The corrections in (1.3)— (1.5) with integer 
powers of n are the analytic corrections to scaling, while A x is the exponent of the leading non-analytic correction 
to scaling. There are additional higher order analytic corrections with exponents —3, —4, • • • , and non-analytic 
corrections with exponents A 2 , A 3 , • • • . The higher order corrections to scaling involve arbitrary mixed terms as 
well, of the general form const. /n J o+-' lAl +J 2A +J 3A3 +'", with non-negative integers ji. In (1.3)— (1.5), AF denotes 
the possible contribution of the so-called anti-ferromagnetic term which occurs for bipartite lattices such as the 
square lattice and the simple cubic lattice. This term alternates in sign, and is sub-leading as the exponent arises 
from polygon-like configurations. 

See [8] for extensive discussions and analysis of the correct asymptotic form for two-dimensional SAW, and [9] 
for a discussion of the full asymptotic form for three-dimensional SAW. 

Due to the importance of SAW, in the absence of an exact solution it is desirable to find simpler models. Recent 
progress has been made in finding interesting solvable models of walks, with non-trivial critical exponents, such 
as prudent walks [10, 11]. Another recent paper [12, 13] takes a different approach, developing a prescription to 
discretize SLE, resulting in a self-avoiding walk which is a finite piece of an infinite SLE curve. 

Other models which are related to SAW are self-avoiding returns and polygons (SAR and SAP respectively) 
and bridges. A self-avoiding return (SAR) is a SAW which is forced to return to the origin, i.e. an n-step SAR 
is identified with an n — 1-step SAW with \uj(n — 1) — uj(0)\ = 1. An SAP is an unrooted unoriented SAR, and 
so there are exactly 2n corresponding SAR for each n-step SAP. Bridges are a subset of SAW with an additional 
condition on the x -coordinates of a walk: uo\ (0) < u)\ (i) < uo\ (n) with < i < n. Both SAP and bridges have the 
same connective constant as SAW, and the same size exponent v as SAW, but different critical exponents for their 
growth in number. 

In this paper we introduce a new variant of SAW, which we call endless self-avoiding walks, or eSAW. We 
will prove that this new model has the same connective constant fx as SAW, and argue that it has the same critical 
exponent v as SAW but with number growth exponent -f c = 1 in any dimension. The number growth has no analytic 
corrections to scaling, and we find strong numerical evidence that the amplitude A c of the number growth is a 
universal quantity. 

We will study eSAW on various two-dimensional and three-dimensional lattices, in particular the square, hon- 
eycomb, triangular, union jack, simple cubic, body centered cubic (bcc), and face centered cubic (fee) lattices. The 
model can be naturally extended to other lattices, and also to other models of walks as described in Sec. 5. 

1.1 Outline 

We will first introduce the endless self-avoiding walk model in Sec. 2, and discuss its expected properties. We 
will then describe the enumeration method in Sec. 3, and analyze the resulting enumerations to estimate critical 
parameters for eSAW in Sec. 4. In Sec. 5 we will discuss various aspects of the eSAW model, and finally conclude 
in Sec. 6. 

2 Endless self-avoiding walks 

The basic idea for constructing the model is extremely simple: we start with the set of self-avoiding walks, and 
endless self-avoiding walks are the subset of walks which are still self-avoiding after they have been concatenated 
with copies of themselves head-to-tail ad infinitum. An example of an eSAW on the square lattice is shown on the 
left in Fig. 2.1, with the root point shown as a solid circle and the concatenated copies as dotted lines, while on the 
right is a SAW which is not an eS AW. Because the walk is extended to infinity, sites on the walk are indistinguishable 
and end-effects are completely eliminated. This definition of eSAW is applicable for any lattice. 



An alternative, illuminating definition of eSAW may be given in terms of words, where we define the walk in 
terms of sequences of steps rather than sites on a lattice. We now specialize to Z d for the sake of concreteness, 
but the argument is similar for other lattices. Let us denote the steps in the positive axis directions as 1, 2, • • • , d, 
and the steps in the negative axis directions as 1, 2, • • • , d. Then SAW of n steps correspond to the set of words 





Figure 2.1: An example eSAW, u ex , on the left; on the right is a SAW which is not an eSAW. 
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which contain no loops, i.e. for which there are no sub-words with equal numbers of 1 and 1, 2 and 2, • • • , d 
and d. Given a walk, the infinite word which results from concatenating that word with itself indefinitely is an 
endless walk. If the infinite word has no loops, then the walk is an endless SAW in correspondence with our earlier 
definition. E.g., if we consider the example SAW from Fig. 2.1, w ox = 121212, then the corresponding infinite word 
is • • • 121212121212121212 • • • . The infinite word has no loops, and hence uj cx is an eSAW. 

When we choose boundary conditions for a lattice model, conventionally this boundary is in real space, i.e. it 
applies to a piece of the lattice. In contrast, we consider eS AW to be SAW with periodic boundary conditions in what 
we could loosely describe as chain space. In chain space, SAW and eSAW are one-dimensional objects, regardless 
of the embedding dimension. 

We are principally interested in the asymptotic behavior of the number of eSAW of n steps, e„, and the square 
end-to-end distance er n (other quantities of interest are defined in Sec. 2.5). In analogy with (1.3)— (1.5) we expect 
that: 

e„~ A^> 7e_1 — + ^l + -" + -^- + "- > ) +AF, (2.1) 
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er n ~ AeA^n 2 ^" 1 [i + J -± + ±± + ... + + ...)+ AF, (2.2) 



e Pn ~ D e n 2v ° (l + - + ^ + -- - + ^ + -- - I +AF, (2.3) 
where ep n = er n /e n . 

In Sections 2.1-2.6 we will specialize to discussing eSAW on 1 d , and will discuss the minor differences arising 
for the honeycomb lattice in Sec. 2.7. 

2.1 Proof that fi e = \i for eSAW on 7L d 

Every bridge of n steps is also an eSAW of n steps, since by construction bridges can be concatenated with them- 
selves (indeed, any bridge can be concatenated with any other bridge and remain self-avoiding). It is also clear that 
every eS AW is in fact a SAW, and hence if B n is the number of bridges of n steps, e n is the number of eS AW of n 
steps, and c„ is the number of SAW of n steps, we must have 

K < e n < c n . (2.4) 
lim BU n < lim e\[ n < lim c^™, (2.5) 
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lim el/ n = M , (2.7) 



i.e. the connective constant for eSAW is precisely the same as that for self-avoiding walks. 
There is strong evidence that b n and c„ obey the scaling laws: 

b n ~ n b -V" (2-8) 
c n ~ V (2-9) 

where b < l 1 . This in turn implies that b < j c < 7; we will argue in the following subsection that in fact -f c = 1. 

2.2 Argument that 7 C = 1 

We will now give a number of heuristic arguments in support of our expectation that 7 e = 1. Our intent is to 
communicate the intuition we have as to why j c = 1. 

'in fact for two-dimensional lattices numerical evidence and SLE arguments provide convincing evidence that b = 9/16. Private communi- 
cation from A. J. Guttmann. 
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Direct renormalization group arguments regarding the partition function for linear polymer models such as SAW 
in two and three dimensions lead to 7 > 1 [15, 16, 17]. Heuristically, this is due to the fact that ends have more 
freedom than the interior of the chain, and the degree of freedom increases as the length of the chain increases. See 
Chapter XL 1.6 of [15] for a qualitative but more detailed argument, and [16] for theoretical arguments regarding 
polymer networks with an arbitrary number of ends. See [17] for a systematic renormalization group calculation 
for linear polymers in conformation space, which explicitly showed the leading order contribution which perturbs 7 
away from 1 is due to an end effect. 

For clarity we reproduce part of the argument of [15] here. Suppose we can express the free energy F of a 
polymer chain in terms of length n and interaction strength u. Then, if we renormalize the chain by a scale factor g 
by integrating out the structure at small length scales, we would expect that we could express the free energy of the 
full chain in terms of the free energy of the renormalized chain as follows: 

F(n, u) = AF(n, u) + F(n/g, m). (2.10) 

We are interested in the behavior of F in the vicinity of the fixed point u* , where u s=s u\ w u* . For SAW, AF 
includes a contribution from interior chain segments, which have two neighbors, and an additional contribution from 
the ends, which only interact with one neighbor. It is the contribution of the ends that perturbs 7 away from 1 . 

For SAW, the enhancement effect becomes negligible for dimensions d > 4 where 7 = 1, while there is a 
logarithmic correction to the number growth in dimension four. For models where the ends have less freedom than 
the interior of the chain, such as bridges and SAP, we would expect (and do indeed observe) that the growth exponent 
is strictly less than unity. 

Now, the complete absence of any end-effects for eSAW should mean that 7 is precisely equal to 1, as there is 
no physical mechanism to perturb the value of 7. For eSAW, AF in (2.10) only includes contributions from interior 
chain segments, and hence j c = 1. 

We may consider eSAW as an object extending to infinity (a rooted self-avoiding chain, in the notation of 
Sec. 2.5). When we directly renormalize a configuration, its global nature remains invariant as it is still extends to 
infinity, and only the local structure of the walk changes. The one large scale alteration is the size of the repeating 
unit, or alternatively the density of boundaries along the walk. For SAW, the sizes of the end-regions of the walks 
scale in a different way to the interior of the chain. That is, for SAW there is an internal length scale (distance from 
the ends) which must be taken into account when the system is renormalized, but for eSAW there is no such length 
scale to worry about. 

The ability to derive a model with 7 = 1 via the imposition of periodic boundary conditions is thus seen as 
fundamentally due to the fact that linear polymers are one-dimensional objects, no matter the dimensionality of the 
space they are embedded in. We therefore believe it is unlikely that other lattice models can be adapted in the same 
way. Of course, this does not preclude the possibility that a new insight could lead to non-trivial two-dimensional 
and three-dimensional models with 7=1. 

We note in passing that the process by which 7 may be altered from its standard value is somewhat subtle: 
theoretical arguments suggested that it should be possible to vary 7 [18] by introducing orientation dependent inter- 
actions into the polymer model. In particular, it was predicted that 7 would be different for SAW on the Manhattan 
lattice [19] than on other lattices due to the fact that SAW cannot be trapped on this lattice (effectively introducing 
an orientation dependent interaction). However, detailed numerical investigation [20] subsequently showed that 7 
was unchanged from its value for SAW on other two-dimensional lattices. 

Another 'endless' walk model is of course the model of self-avoiding returns. The number of self-avoiding 
returns grows as n a ~ 2 /j, n , and the hyperscaling relation dv = 2 — a can be interpreted in the following way. SAW 
with free end-points have typical size n" ', and so the fraction which return to the origin should be 0(n~ dv ). Thus 
we should expect the number of SAR to be 

#(SAR) ~ n a n- du n n , (2.11) 

where wT dv is the factor due to geometric restriction, and n a is an enhancement factor which parallels the n 7o_1 
factor for eSAW. There is strong theoretical and numerical evidence in support of the hyperscaling relation, which 
implies that a = 0. Since there are no geometric restrictions for eSAW, we regard the validity of the hyperscaling 
relation for self-avoiding returns to provide support for our argument that 7 e = 1. 

We readily acknowledge that our arguments for 7 e = 1 fall far short of a proof, and will verify 7 e = 1 to good 
precision via series analysis methods in Sec. 4. 

2.3 Argument that u e — v 

eSAW are different from SAW in that they have no ends, and so in a small region of the chain near the end one may 
expect that, on average, eSAW and SAW are qualitatively different. However, the size of SAW and similar objects 
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is dominated by the region of the chain far from the end, and this is why SAW, SAP, bridges, and eSAW all have the 
same size exponent v. 

More explicitly, following Caracciolo et al. [21], we define 

£ = ^Pn, (2.12) 

e =£L (2.13) 

Sn 

Then in the limit n — > oo, |r| — » oo, g fixed, we expect that the end-to-end distribution function for SAW is [22, 21] 

P{g) = ^f{g){l + 0{n-^)) (2.14) 

where for large g 

f(g)~ focQ a c W (-Dg s ), (2.15) 

with a and S both positive. Now, eSAW are a subset of SAW, and so if we suppose that v e > v, then (2.15) would 
imply that e„ could not grow asymptotically as e„ <~ n a [i n . The most reasonable interpretation, based on the 
assumption that e n <~ n a /i™ and ep n <~ n 2iy °, is that f e < 

In addition, we expect that eSAW which return close to the origin are suppressed, i.e. polygon-like configurations 
are relatively rarer amongst eSAW than they are amongst SAW. The endless condition induces an effective repulsion 
between the end-points, and thus we would expect that v c > v. Therefore we conclude that v c = v. 

As an aside, we expect that the end-point repulsion would result in a larger amplitude for the mean square 
end-to-end distance for eSAW than for SAW. 



2.4 Analytic corrections to scaling 

For SAP on two-dimensional lattices, there is strong numerical evidence that the non-analytic corrections to scaling 
fold into the analytic background term [23, 24, 8]. 

For eSAW, we expect that the analytic corrections to scaling will fold into the analytic background term, due to 
the fact that 7 e = 1. If we define the generating function for the number of eSAW as 

oo 

C (x) = ]Te„x", (2.16) 

n=l 

then we expect that the asymptotic behavior of the e„ series is determined by the singularity or singularities of 
C e (x) which are closest to the origin. Neglecting the sub-leading anti -ferromagnetic singularity (which may have 
zero amplitude for eSAW in any case; see Sec. 5.3), we expect that 

C e (x) = A{x) + B(x)(l - x/x c )-* (1 + c(l - x/x c ) Al +•••), (2.17) 

where A{x) and B{x) are analytic functions. If we write 

B(x) =A e + p 1 (l- x/x c ) + (3 2 (1 - x/x c f + /3 3 (1 - x/x c f + ■■■ , (2.18) 

then, under the assumption that 7 = 1, all of the higher order terms with coefficients 02, fiz, • • • will be absorbed 
into the analytic background term A(x)l 

Thus, in (2.1) we expect that a\ = a 2 = • • • = 0. We will test this hypothesis, as there could conceivably be a 
confluent logarithmic correction which would result in non-zero coefficients aj. 

2.5 Irreducible eSAW and self-avoiding chains 

The more natural definition of an endless SAW (to us), which we will call a self-avoiding chain (SAC), is related to 
eSAW in the same way that self-avoiding polygons are related to self-avoiding returns. SAC are equivalence classes 
of eSAW which correspond to the same 'shape', regardless of the root point, or the direction of traversal. 

For an n-step eSAW on 7L d there are n possible root points, and two directions of traversal, suggesting that the 
number of eSAW of n steps should be 2n times the number of SAC. However, there is an additional complication for 
the eSAW / SAC correspondence which does not occur for self-avoiding returns / self-avoiding polygons: eS AW can 
have repeating sub-units which makes the infinite chain corresponding to a given eSAW non-unique. We illustrate 
this in Fig. 2.2, where it can be seen that the same infinite chain results for the two eSAW shown. We call those 
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eSAW with no repeating units irreducible eSAW. The number of eSAW, e„, and the square end-to-end distance r n 
can then be straightforwardly expressed in terms of the number irreducible eSAW, i n , and the square end-to-end 
distance of irreducible eSAW, ir n : 



,2 



En 
*2 



2 lr i- 



The inverse expressions can be obtained, if desired, via a Mtibius transform involving the Mobius function /i: 



in = Yl I 



E( n\ n 2 
Hi J .r" V 



(2.19) 
(2.20) 

(2.21) 
(2.22) 



An alternative interpretation is that the word corresponding to an irreducible eSAW cannot be decomposed into 
identical, repeating sub-words. 




Figure 2.2: Examples of two eS AW which correspond to the same infinite chain. The eSAW on the left is irreducible. 

Note that one does not need to worry about "irreducible" SAR configurations precisely because of the self- 
avoidance condition. If we were to consider weakly as opposed to strictly avoiding self-avoiding returns / polygons 
instead then the same complication would occur, as loops could then repeat. 

Now for each irreducible eS AW of n steps there are n possible root points, as can be seen for our example eSAW 
in Fig. 2.3. In addition there are two possible orientations for the walk, and hence there 2n irreducible eSAW for 
each self-avoiding chain. If we denote the number of SAC with an irreducible unit of n steps as d n , then 

*„ = 2nd n . (2.23) 

In terms of words, converting from rooted to unrooted eSAW can be done by enforcing an ordering on the set of 
words corresponding to a given chain, and then choosing a single representative of this set. A natural way to do this 
is to use lexicographic ordering by placing an order on the steps themselves, e.g. 1 < 1 < 2 < 2. In Fig. 2.3 this 
results in the ordering 121212 < 121212 < 121212 < 212121 < 212121 < 212121; the eSAW in the top left of the 
figure, 121212, would thence be chosen as the representative. Orientation independence on words can be enforced 
by defining words which are mapped to each other under the involution IhI and 2 <H- 2 to be equivalent. 

In summary, endless walks may be reducible or irreducible, rooted or unrooted, and oriented or unoriented. 
Using this taxonomy, the objects we will study in this paper are shown in Table 1. We expect that the asymptotic 
behavior of the irreducible and reducible series will be exactly the same, i.e. i n ~ e n and ir n ~ er n . Since 
d n = *n/(2n), the exponent for the growth in d n is (7 e — 2) instead of (7 e — 1), but otherwise the asymptotic 
growth for d n is given in (2.1). 

We show the self-avoiding chains with irreducible units of 1-4 steps in Fig. 2.4. 

2.6 Properties of eSAW 

We will now briefly mention some properties of eSAW which may not be immediately apparent. 

Firstly, although all bridges are eSAW, the converse is not true. For example, the eSAW in Fig. 2.5 is not a bridge, 
and no choice of start and end points on the self-avoiding chain corresponding to the eSAW results in a bridge. A 
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Figure 2.3: Irreducible eSAW which are equivalent up to the root point. 



Object 


Properties 


Sequence 




reducible, rooted, oriented 




eSAW 


square end-to-end distance of eSAW 


er n 




mean square end-to-end distance of eSAW 


ep n 




irreducible, rooted, oriented 




irreducible eSAW 


square end-to-end distance 


ir n 




mean square end-to-end distance 


ipn 


self-avoiding chains 


irreducible, unrooted, unoriented 


d„ 



Table 1: Definitions of objects and their corresponding sequences. 



necessary and sufficient condition for a SAC to allow for the choice of a bridge as the repeating unit, or equivalently 
for the word representation of an eSAW to have a cyclic permutation which is a bridge, is that there must exist a 
vertical line which cuts only a single bond of the SAC. The proof of this statement is obvious by inspection. 

Going further, a natural generalization of bridges would be to allow for directions other than vertical for the 
cutting line. If all SAC were to admit at least one straight line which cuts only a single bond, then we would regard 
eSAW as a generalization of bridges. However, it is straightforward to derive eSAW for which any straight line must 
past through at least two bonds of the SAC, for example by slightly perturbing the eSAW in Fig. 2.5. 

Trapped self-avoiding walks are SAW which cannot be extended indefinitely from (at least) one of the ends and 
still remain self-avoiding (see e.g. [25]). By definition, no trapped self-avoiding walk is an eSAW. Conversely, every 
non-trapped self-avoiding walk is a sub-walk of an infinite number of eSAW, since the walk can always be extended 
far away from any problematic regions. Equivalently, all eSAW are Kesten patterns [26], and we expect that all 
patterns have positive density in sufficiently long eSAW. Examples of trapped and non-trapped SAW are shown in 
Fig. 2.6. 

Unlike SAW, eSAW in dimensions d > 2 allow for true knots to occur when we consider them as infinite chains. 
An example knotted eSAW is shown in Fig. 2.7; it corresponds to the word 11233221333223321, and is of length 
17. This compares to the shortest knotted SAP on the simple cubic lattice, which has 24 steps [27, 28]. We do not 
know if 17 is the minimum length for an eSAW knot on the simple cubic lattice. 

Finally we note that in two-dimensions all eSAW must have zero winding angle. The winding angle has featured 
in recent proofs of exact results for self-avoiding walks [7, 29]. 

2.7 eSAW on the honeycomb lattice 

The set of direction vectors for the neighboring sites of the honeycomb lattice depends upon which of two bipartite 
sub-lattices one is currently on. This contrasts with the other lattices considered in this paper which possess full 
translation invariance. As a consequence, our definition of eSAW implies that on the honeycomb lattice eSAW may 
only have an even number of steps. An example of an eSAW for the honeycomb lattice is shown in Fig. 2.8. This 
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Figure 2.4: The self-avoiding chains with irreducible units of 1-4 steps on the square lattice. By inspection, d\ = 2, 
d 2 = 2, d 3 = 4, and c? 4 = 8, in agreement with the values in Table 2 of A. 




Figure 2.5: An eSAW for which no cyclic permutation is a bridge. 

in turn means that a 2n step self-avoiding chain on the honeycomb lattice has n possible root points, rather than 2n. 
The relations between the different sequences must be modified accordingly: 

e 2n = J2 i2 i> (2 ' 24) 
n 2 



er 2 n = ^2^ir2j, (2.25) 

i2n = ^2n(-)e 23 , (2.26) 

ir2n = ^2^(^) r ^er 2j , (2.27) 

i 2n = 2nd 2n , (2.28) 



where once again /j, is the Mobius function. 
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Figure 2.6: In solid lines, a trapped SAW on the left and a non-trapped SAW on the right (a possible continuation is 
shown as a dotted line). Neither of these are eSAW, but the non-trapped SAW may be a sub-walk of an eSAW. 



I rrf'TlJI fu 
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Figure 2.7: A knotted eSAW. 



It is possible to broaden the definition of eSAW to allow for eSAW with an odd number of steps by altering 
the boundary conditions, e.g. by interchanging left turns and right turns at the concatenation point. However, this 
introduces some technical difficulties 2 , and in the absence of a pressing reason to do so we will not explore this 
option further here. 

One key advantage for the analysis of eSAW on the honeycomb lattice is that the connective constant is known 
to be exactly [5, 7]. We will exploit this fact later in our analysis in Sec. 4. 



/ **-* 




Figure 2.8: Example eSAW on the honeycomb lattice. Note that the repeating unit has an even number of steps; if 
one step were removed then the eSAW would no longer be able to remain on the honeycomb lattice. 

The arguments for v c and -f c are not altered for the honeycomb lattice. However, the proof that [i c = [i would 
need to be modified. We will not prove that [i c = [i here, but are confident that this statement is indeed correct. 



3 Enumeration method 

We performed enumerations on the square, honeycomb, union jack, simple cubic, bcc, and fee lattices using a basic 
backtracking algorithm. Backtracking works by explicitly generating all of the combinatorial objects of interest by 
recursively building up these objects from smaller objects. 

We will not describe backtracking in detail, but note that it has been used to study SAW and related models 
for over 50 years (see section 7.3 of [30] for references). Here we focus on the differences between our algorithm 
for enumerating eSAW and the standard backtracking approach for enumerating SAW. More details on how to 
implement a backtracking algorithm for the enumeration of SAW are given in [9]. 

2 E.g., in this context it seems more natural to use 'mid edges' as in [7] to define an eSAW with twisted boundary conditions, because a single 
edge does not carry any information about whether a left or right turn comes next. 
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The key difference between backtracking enumeration of eSAW and SAW, is that as an eSAW is built up step by 
step one has to check for self-intersections for the whole infinite chain. This is apparently quite a non-local effect, 
but it can in fact be converted to a local interaction provided the end-to-end vector is known. 

Suppose we have end-to-end vector v, then if we know that x is a site on the walk then necessarily all sites 
y = x + jv, j G Z must also be part of the endless walk. We should therefore test for self-intersections for an 
endless walk uj by checking for self-intersections in reduced coordinates cj(i)/v. 

Our algorithm generated all eSAW with a given end-to-end vector, v, and used reduced coordinates to test for 
self-intersections. Starting at the origin, we incrementally grow our walk by considering all possible ways of adding 
an additional step which are compatible with the end-to-end vector. As the walk is extended we add uj(i) to the list 
of visited sites, and we check to see ifw(«)/v = uj(j)/w for any j < i. If there is an intersection we know that this 
sub-walk can never be a valid eSAW, and so we truncate the backtracking tree. 

The reduced coordinates of visited sites were stored in a hash table to ensure that intersection testing took 0(1) 
CPU time. Optimal backtracking implementations take CPU time which is of exactly the same order as the number 
of objects being generated, and in this sense our implementation was optimal. It certainly would be possible to 
obtain a small constant factor improvement in performance through fine tuning the implementation. 

We performed our enumerations on the following lattices in two dimensions: square, triangular, honeycomb, 
and union jack, and in three dimensions: simple cubic, body centered cubic (bcc) and face centered cubic (fee). The 
two-dimensional lattices are shown in Fig. 3.1. The lattice we refer to as "union jack" is non-standard, as it allows 
for crossings, and may be properly described as a square lattice with additional bonds connecting next-to-nearest 
neighbors [31, 32]. 




Figure 3.1: Two-dimensional lattices, from left to right: square, triangular, honeycomb, and union jack. 

We enumerated e„ and er n directly, and the fact that the resulting value for d n must be an integer places quite 
a strong consistency condition on our enumerations. Another condition that our enumerations should satisfy comes 
from the observation that eSAW form equivalence classes under transformations of the underlying lattice. For 
the square lattice, eSAW which explore both available dimensions must form an equivalence class under the full 
symmetry group, 0(2), which means that e„ — 4 must be divisible by 8. The corresponding condition for the simple 
cubic lattice is that e$f mple cubic) - 3e^ qimre) + 6 must be divisible by 48, the number of elements of 0(3). We have 
confirmed that each of these conditions hold for our data. 

Interestingly, SAC do not necessarily exhibit the full symmetry group of the underlying lattice. The symmetries 
under which SAC are invariant is determined purely by the end-to-end vector. 

The combination of integer values for d n , correct symmetry factors for e„, and correct low order values gives us 
confidence that the enumerations reported in Tables 2-4 in A are correct. As an aside, we note that for n prime, the 
sequences must satisfy 



We performed our enumerations on a SunFire X4600M2 with 2.3GHz AMD Opteron CPUs. For each lattice we 
divided the backtracking enumeration amongst 20 different processors. The total number of CPU hours required for 
the enumerations was: 12000 for the square lattice, 2100 for the triangular lattice, 2400 for the honeycomb lattice, 
2000 for the union jack lattice, 5600 for the simple cubic lattice, 7000 for the bcc lattice, and 1300 for the fee lattice. 

Our enumeration data are collected in Tables 2-4 in A. Note that we only report r„ for the square and simple 
cubic lattices. A bug in our enumeration program for the other lattices meant that the data we collected suffered from 
rounding errors. These data would allow us to estimate the amplitude D c , but since it is a non-universal quantity we 
regard D e to be of relatively minor interest. We will present enumerations of r n for other lattices in future papers. 

Finally, note that we do not regard basic backtracking as best practice for performing the enumerations, rather 
as a sensible first step towards studying this new model. We briefly mention alternative enumeration approaches in 
the following subsection. 



C-2n = l2n + «2 = 2nrf 2 n + 2d 2 

e„ = i n + h = 1nd n + 2dx 



honeycomb lattice, 
all other lattices. 



(3.1) 
(3.2) 
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3.1 Alternative enumeration methods 



More sophisticated enumeration algorithms hold promise for future progress, but implementing them for eSAW will 
require overcoming certain challenges due to the apparent non-locality of eSAW self-interactions. Our intuition is 
that this may present a significant technical challenge for the approaches mentioned below, but not an insurmountable 
obstacle. 

The most straightforward approach to improving our enumerations would be to implement the two-step method [9] . 
This could certainly be used to extend the series in A by a few terms in each case for the same computational effort. 
We are not sure if the lace expansion technique described in the same paper could be useful for enumerations of 
eSAW. 

The most powerful known method for three-dimensional walks is the length-doubling algorithm [33, 34]. It is 
not clear to us how this exciting new enumeration algorithm could be applied to eSAW. 

We believe the most promising approach is the finite lattice method [35, 36], which has been applied to the 
problem of enumerating SAW [37] and SAP [38] on two-dimensional lattices. A recent improvement to the finite 
lattice method [38] raises the possibility that the method could also be usefully applied in three dimensions. 

4 Analysis 

Our goal in analyzing the series from Appendix A is to extract as much information as possible. In particular we 
wish to provide numerical support for our argument that 7 = 1 in Sec. 4.3 and u c = v in Sec. 4.4. In Sec. 4.5 
we provide evidence for the surprising fact that the amplitude A c is universal, while confirming that there are no 
analytic corrections to scaling for the number growth. 

For our analyses we use arguably the two most powerful, general purpose methods of series analysis available: 
the methods of direct fitting and differential approximants. We succinctly describe each of these methods below in 
Sections 4.1 and 4.2 respectively. 

In our analyses we will ruthlessly exploit our prior knowledge of critical parameters, in particular the values 
of the connective constants of the various lattices. For two-dimensional lattices, the best estimates for \i come 
from series derived from the finite lattice method. For the square lattice, /j, = 2.638158 530 35(2) [38], while 
for the triangular lattice /j, = 4.150 797226(26) [39]. For the honeycomb lattice we have the exact result fi = 
\Jl + \71 = 1.847 759 065 022 • • • [5, 7]. For the union jack lattice there are no accurate prior estimates of the 
connective constant, and so in contrast to the other series we perform unbiased fits. From a recent Monte Carlo 
calculation, we know that \i = 4.684 039 931(27) [40] for the simple cubic lattice. Using the length-doubling 
algorithm, Schram [41] has derived series for the bcc and fee lattices to 28 and 24 steps respectively. A preliminary 
analysis of these data leads to fj, = 6.530 521(6) (bcc) and fj, = 10.037067(7) (fee). 

The values above for the connective constants are in all cases sufficiently accurate that we can consider them 
to be "exact", i.e. we do not have to measure the effect of the error in the estimates of \i on our fits. This is not 
usually the case with biased fits. Our situation is special because this knowledge comes from other sources, and \i is 
known to much higher precision than could possibly be obtained from our short eSAW series. We have tested this 
assumption for ourselves, but will not discuss this issue any further in the analysis below. 

We found that eSAW series are in every case somewhat better behaved than the irreducible eS AW series. Gener- 
ically, they are quite different to begin with but then the fits converge for larger n, since they are asymptotically 
identical. We attribute this to the fact that for small n the effect of repeating sub-walks for composite n values is 
quite high, whereas when n is large this effect becomes negligible. To reduce visual noise only the eSAW series are 
plotted. 

We are mindful of the following observation due to Guttmann [42] with reference to extrapolations from series 
analysis: "error bounds are generally referred to as (subjective) confidence limits, and as such frequently measure 
the enthusiasm of the author rather than the quality of the data.". Where we give an explicit confidence interval for 
our estimates we will try to moderate our enthusiasm. We will also present all of our fits in graphical form so that 
readers can make their own judgments. 

4.1 Direct fitting 

The method of direct fitting involves fitting the asymptotic form for the series with the series coefficients. This 
approach is especially useful for finding amplitudes. 

The idea is to truncate (2.1)-(2.3) at some finite order, and then fit a subset of coefficients from the corresponding 
series to this truncation. Rather than going in to any great depth about the general method, we will instead explain 
exactly what fits are performed in the analysis below. The most detailed explanation is in Sec. 4.3. 

See [42, 9, 43] for detailed explanations of the direct fitting approach. 
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4.2 Differential approximants 



In the method of differential approximants, the unknown generating function is represented by the solution to an 
ordinary differential equation of the form: 



' d 



i=0 



The functions P and Qi are polynomials, of degrees L and N. We choose < L < 10, Nk > 3 (which guarantees 
at least three regular singularities), and take Qk to have highest-order coefficient equal to 1. We choose K — 2 for 
the two-dimensional lattices, and K = 1 for the three-dimensional lattices as these series are generally shorter. The 
order of the polynomials was chosen so that \N — Nj\ < 3. Given coefficients do, . . . , a^, the polynomials P, Qi 
are chosen so that the polynomial J2n=o a « z ™ s °l ves the differential equation to within an error of order z N+1 . This 
choice is made by solving a system of linear equations in L + K + 1 + X^i=o ^* unknowns, determined from N + 1 
known coefficients. 

Once the approximant for a given choice of L, K, and Ni has been determined, one may then find the singular- 
ities of the differential equation by determining the zeroes of Qk{z). The exponent corresponding to each regular 
singularity is then determined from the indicial equation. For the singularity z i7 the corresponding exponent is: 

K = K-l-^^. (4.2) 
ZiQ'{zi) 

The physically relevant singularity is the nearest singularity on the positive real axis. If we denote this singularity 
as z c , then the connective constant fi = l/z c . We denote approximants which have singularities near the physical 
singularity, or clearly incorrect singularities on or near the positive real axis, as "defective". To estimate z c and the 
corresponding critical exponent we determine the leading singularity for a large number of choices of L and Ni, 
and prune defective approximants away. We then plot these approximant estimates of the critical exponent against 
n (the power of the highest order series coefficient used) in order to extrapolate and obtain a final best estimate for 
the critical exponents 7 and v e . 

If desired, it is straightforward to ensure that there is a singularity at a biased value of z c by introducing an 
additional linear equation: 

Qk(zc) = 0. (4.3) 

In this paper, we are able utilize accurate estimates for /x from longer SAW series, and so all of our differential 
approximant analyses are biased. 

See Guttmann [42] for far more information about differential approximants and series analysis, or Guttmann 
and Jensen [43] for a recent overview. 



4.3 7e 

To determine 7 C for the square, honeycomb, and triangular lattices we take the logarithm of e n /fi n , and note that 
the asymptotic form we derive from (2.1) is 



log (e„//i n ) ~ (7c - 1) logn + log A c + 



n n Al 



(4.4) 



We single out the leading order neglected terms by placing them in large square brackets. Note that we have 
substituted [i for fi c due to the proof in Sec. 2.1. We then solve for 7 e and log A e by using terms from the sequence. 
For the triangular and honeycomb lattices we use consecutive terms, say ek-i and for the triangular lattice and 
e fe _ 2 and e k for the honeycomb lattice. For the square lattice we find that we obtain smoother fits by fitting alternate 
terms e^-i and e k - We label the resulting estimates for 7 e as 7 e ,fc, and plot these estimates in Fig. 4.1. Note that 
we choose to take the logarithm in the above equation simply as a matter of convenience, so that the equation to be 
solved for 7 and log A c is linear. 

Now, we expect that the error term in our 7 C „ fits is of the order of the first neglected term. The asymptotic 
behavior of two-dimensional SAW is well understood [8], thanks to a combination of series analysis of long series 
generated via the finite lattice method, and high powered Monte Carlo simulations. This leads us to believe that 
Ai = 3/2 exactly. The argument of Sec. 2.4 suggests that ai = 0, and so we expect that plots of 7 C „ versus n~ Al 
should be linear, and this is confirmed for each lattice in Fig. 4. 1 . We confirmed that plots of 7 6i „ versus 1 jn possess 
significant curvature, indicating that a\ must be close to zero in each case; we will estimate a\ in Sec. 4.5. 



12 



1.000 
0.995 
0.990 







- * ' ' x ' 

x-x' x "' x " 
. xx-x* x * 


X 

- - - ' ' X 


1.01 
























o- - -O. - 


1.00 

7e,n 










*-* 




0.99 






X 


honeycomb 




0.98 






O 


square 


-~.A . 






O simple cubic 


1 1 1 


~~ A 

* 


triangular 
union jack 
1 1 


* 

II 


A- 

1 


0.97 


• bcc 

" A fee 

1 1 



0.000 0.004 0.008 0.012 0.00 0.05 0.10 0.15 0.20 0.25 

(nlog/i)" 1 - 5 (rclog/i)" - 528 



Figure 4.1: Direct fit estimates of 7 for two- and three-dimensional lattices (left and right respectively). 



We employ one additional trick in our plots of direct fit estimates. For a lattice with connective constant /i, when 
we plot against a remainder of the form rT a we normalize it as x — {n\og[i)~ a . This ensures that for any given 
value of x the number of eSAW is approximately the same, independent of the lattice. Loosely speaking, we want 
the x coordinate to be correlated with the information content of the series. 

We also obtain direct fits for 7 C from the union jack lattice, solving 

log e„ ~ n log pi, + (7e - 1) log n + log A c + 

for log [i, 7c, and log A c using coefficients efc_ 2 , ek-i, and e/j. These estimates are not as well behaved because we 
do not have the benefit of a prior high precision estimate for /x. 

Examining Fig. 4.1 it is clear that all of the estimates for 7 o n are converging smoothly to 1, as predicted. As 
a visual guide we fit a straight line to the five highest order estimates for each lattice, and their y-intercepts are 
extremely close to 1, especially for our high precision data from the square, honeycomb and triangular lattices. 
From the figure, j c must be within the range 1.000(1). We could do better with the direct fits, but we will see that 
the differential approximant analysis below is more than adequate for the task of placing a tight bound on -f c . 

For the three-dimensional lattices, we once again fit the asymptotic form 

log (e n /fi n ) ~ (7c - 1) logn + log A c + 

For three-dimensional SAW, the asymptotic form is well understood theoretically, but the numerics are far more 
difficult. In [9] there is extensive discussion of the asymptotic form, while in [44] there is the most detailed and 
accurate numerical investigation of corrections to scaling that is currently available. For three-dimensional lattices 
we have Ai = 0.528(12) [44] as the dominant neglected term, but in other respects we proceed in the same manner 
as for the two-dimensional lattices. For the simple cubic and bcc lattices we fit ek-2 and e&, while for the fee lattice 
we fit efc^i and e/j. 

Our fits for 7 C are shown in Fig. 4.1. It is clear in each case that the estimates are smoothly converging to 7 = 1. 
The fits as shown suggest that 7 = 1.005(10), which encompasses 7 = 1 as predicted. We do need to consider the 
effect of the known error in Ai = 0.528(12), and so re-performed the analysis for Ai = 0.516 and Ai = 0.54. In 
our case, because Ai is known to high precision from another source, this does not materially change our estimate 
for 7c 

As a technical aside, we find that for a number of lattices, particularly the square lattice, the simple cubic lattice, 
the bcc lattice, and to a lesser extent the honeycomb lattice, there is quite a high degree of oscillation with our direct 
fit estimates. We reduced the amplitude of the oscillations by fitting alternate terms, and on some occasions we also 
smoothed the fits by averaging consecutive estimates. 

As mentioned in the description of the asymptotic forms (1.3)— (1.5), for SAW there is a contribution to the 
asymptotic form from the so-called anti-ferromagnetic singularity. In the analyses performed in this paper we avoid 
having to deal with this singularity directly by treating the series arising from odd and even terms separately. As 
explained in Sec. 5.3, there is no clear evidence as to whether the oscillations in the series are strictly odd-even, 
which would be a signal of the anti-ferromagnetic singularity on the negative real axis. The approximate alternation 
of some series is not conclusive, as it may be due to singularities near to but not on the negative real axis. Instead, 
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(4.5) 



n n Al 



(4.6) 
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our choice to average consecutive terms and fit alternating terms is based purely on the observation that this makes 
the fits easier to extrapolate. 

We now proceed to estimate 7 C for the various lattices via the differential approximant method described in 
Sec. 4.2. Our fits are reproduced in Figures 4.2, 4.3, and 4.4. 
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Figure 4.2: Differential approximant estimates of 7 on the square (left) and triangular (right) lattices. 
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Figure 4.3: Differential approximant estimates of j e on the honeycomb (left) and simple cubic (right) lattices. 
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Figure 4.4: Differential approximant estimates of j c on the bcc (left) and fee (right) lattices. 
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For the two-dimensional lattices, the differential approximants appeared to converge better for the square and 
honeycomb lattices. Nonetheless, combining information from all three suggests that 7 = 1.000(1) in line with 
our direct fit estimate. 

For the three-dimensional lattices the estimates are not as well converged, but this is to be expected because the 
series are shorter, and the leading correction to scaling is larger. Nonetheless, the trend in each case seems to be 
converging to one, and suggests that 7 C = 0.99(2). We acknowledge that extrapolating such estimates is a fraught 
process, especially since there is clearly still systematic drift in the estimates. In addition, we already 'know' the 
answer, and this can easily introduce a subconscious bias. However, our interest here is less in deriving an estimate 
for 7e, than in showing that 7 C = 1 is the most plausible interpretation of the available evidence. 

We consider the series analyses presented to be quite convincing in their support of our prediction that 7 e = 1 for 
eSAW on two- and three-dimensional lattices. For the reminder of the analysis section we will assume that -f c = 1, 
and use this information to bias our fits. 



4.4 u e 

For the purpose of estimating u c we find that the er„ series is far more useful than the ep n series. We only have 
reliable data for the square and simple cubic lattices, but we given that v e is a universal quantity this is not a problem. 

We assume that 7 C = 1, and substitute this into the asymptotic expression for er„, and then perform direct fits 
on the following asymptotic form: 



log(er„/ M n ) - {2v e + 7e - 1) logn + log(A c L> c ) 
- 2^ logn + log(A £> c ) + 



/i 



So 



n n 



Ai 



h _9o_ 
n n Al 



(4.7) 
(4.8) 



For each of the lattices we fit alternating terms, ek-2 and efe, to obtain estimates for v e ^ and the amplitude, and we 
plot the results in Fig. 4.5. 
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Figure 4.5: Direct fit estimates of v c for the square (left) and simple cubic (right) lattices. 



For the square lattice, it is very clear that the leading analytic correction to scaling has non-zero amplitude, as 
shown by the linear behavior of the fits in Fig. 4.5. The convergence to v c = v = 3/4 is extremely convincing; we 
conservatively interpret our results as suggesting v c = 0.7500(1). 

For the simple cubic lattice, the convergence is not as good, but this is to be expected for such a short series 
with strong corrections to scaling. Our estimates in Fig. 4.5 suggest that v c is in the vicinity of 0.586, while the best 
estimate for v, v = 0.587 597(7) [44], is shown on the plot. 

The results from our differential approximant analyses are shown in Fig. 4.6. Technically, we are in fact plotting 
estimates of v c + (j c — l)/2, but we assume that 7 = 1 as stated previously. 

For the square lattice the convergence of the differential approximants is very good, and v c = 0.7500(1) is a 
conservative interpretation of the data. 

For the simple cubic lattice, once again there is still a systematic drift occurring at n = 24, and the estimates for 
v c are approximately clustered around 0.590-0.593. This is consistent with v c = v, but not completely convincing. 
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Figure 4.6: Differential approximant estimates of u c on the square (left) and simple cubic (right) lattices. 

We believe that the numerical evidence supporting v c = v is unequivocal for the square lattice. For the simple 
cubic lattice we cannot be so definite due to the short series available, but we nonetheless strongly believe that 
v c = v is the most plausible interpretation of the evidence for three-dimensional lattices also. 

To obtain definitive evidence in support of v c = v in three dimensions we will either need to significantly extend 
the series, or apply the pivot algorithm [44]. The pivot algorithm in the context of eSAW is discussed in Sec. 5. 



4.5 A c 

Assuming that j c = 1, the asymptotic form for e„ is given by: 



i + — + -4 + 
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a i 
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n 
a 2 



n A 2 



n A 2 



For the square, honeycomb, and triangular lattices we perform "direct fits" to 

~A c ai A e b ~ 



en/f/ 1 



A e + 



In this case the fits are trivial: one merely has to calculate the ration e n //i™. For the union jack lattice, we fit 



log e„ <~ n log (i + log A c + 



ai b 
n n Al 



(4.9) 
(4.10) 

(4.11) 

(4.12) 



We plot our estimates against l/n Al in Fig. 4.7, noting that once again this gives clear linear trend which suggests 
that the analytic correction has zero amplitude. 

We now observe a remarkable fact: our amplitude estimates appear to be converging towards the same constant. 
The coincidental convergence of the estimates from series from four different lattices seems highly unlikely, and 
thus we conclude that A c is a universal quantity! 

We test this conjecture further by assuming that a = 0, and fitting 



A c + ^ + 

n Zii 



A e a 2 
n 2 



(4.13) 



for the square, honeycomb, and triangular lattices only. We plot the resulting fits for A c against the assumed form 
of the next-to-leading non-analytic correction, i.e. 1/n 5 / 2 [8]. Although we expect that a 2 = we do not rely upon 
this fact; essentially the same estimates are obtained for A c from plots against 1/n 2 . 

The resulting fits show a high degree of convergence. The scatter between the different lattices is consistent with 
there being non-trivial sub-leading corrections to scaling, which is not surprising for such short series. 

We estimate that A c = 1.57075(10) for two-dimensional eSAW. We note in passing that tt/2 = 1.570796326794- 
lies within the confidence interval of our estimate, although we have no theoretical basis for conjecturing that 
A c = tt/2. 
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Figure 4.7: Direct fit estimates of A e for two-dimensional lattices. On the left is the result of fitting asymptotic form 
with no correction to scaling, while on the right the leading correction to scaling is fitted. 



For the simple cubic, bcc, and fee lattices we perform direct fits to 
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(4.14) 



(4.15) 



For the simple cubic and bcc lattice we once again fit alternate values, while for the fee lattice we fit consecutive 
values. We plot the resulting fits for A c in Fig. 4.8. We plot the fits from (4.14) against the leading neglected term, 
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Figure 4.8: Direct fit estimates of A c for three-dimensional lattices. On the left is the result of fitting asymptotic 
form with no correction to scaling, while on the right the leading correction to scaling is fitted. 

l/n Al . For the fits from (4.15), even if the analytic correction to scaling in has zero amplitude there are still two 
competing non-analytic corrections with exponent in the vicinity of —1. Thus we plot fits from (4.15) against 1/n. 

Although we do not show these fits, we have tested the effect of using values of Ai = 0.528(12) within its error 
bar. In the fits from (4.15) this has the effect of extending the extrapolations to the vicinity of 1.180 and 1.185. 

The estimates for A c from all three lattices once again appear to be converging to a single value to good precision. 
This strongly suggests that A c is a universal quantity for three-dimensional eSAW also. We estimate that A c = 
1.183(3) 

Finally, we will attempt to measure the amplitude of the leading analytic correction to scaling for two-dimensional 
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eSAW, which we expect to be identically zero. To do so, we fit 



e n /,^A e+ A ^ + A ^ + 
n n^ 1 



A e a 2 A e bi 



(4.16) 



for the square, honeycomb, and triangular lattices only. The inclusion of the Ai non-analytic correction is necessary 
to be able to obtain sensible estimates for a\. We plot the resulting estimates for a\ from all three lattices in Fig. 4.9. 
We plot the fits against l/n A2_1 since this is the relative size of the first unfitted term (assuming that a 2 = 0). 
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Figure 4.9: Direct fit estimates of the amplitude of the leading analytic correction to scaling, a\, for the two- 
dimensional lattices. 

We can see that for all three lattices the estimates for a\ are in the vicinity of zero, in agreement with our earlier 
observations based on the trends of our fits. The honeycomb lattice fits are still quite scattered, but the fits for the 
triangular and square lattices each suggest that a\ = 0.00(2). Although this is perhaps not as strong numerical 
evidence as that in favor of the universality of A c , we nonetheless think that the evidence is convincing. 

We do not attempt a similar analysis for the three-dimensional lattices, because there is no realistic possibility 
of obtaining meaningful fits for a\. We discuss this further in Sec. 5.2. 



5 Discussion 

In the previous section, we confirmed that indeed j e = 1 and v c = v as we had predicted. 

We now discuss our discovery that A c appears to be universal in Sec. 5.1, the nature of corrections to scaling for 
e n in Sec. 5.2, and the nature of the anti-ferromagnetic singularity for eSAW in Sec. 5.3. Then we mention some 
open questions for eSAW which will require further study in Sec. 5.4, before discussing some possible applications 
ofeSAWinSec. 5.5. 



5.1 Universal amplitude 

In Sec. 4.5 we presented strong numerical evidence that the amplitude for the number growth of eSAW is a universal 
quantity, with A c = 1.57075(10) in two dimensions, and A c = 1.183(3) in three dimensions. We had not antici- 
pated this a priori, unlike the case for our argument that j e = 1. Indeed, the eSAW model was constructed with the 
goal of creating a model with -f c — 1. We do have a heuristic post hoc argument to explain why the amplitude is 
universal, which we present below. 

Once ends have been eliminated, p, is simply a measure of local flexibility, and the polymer then only has one 
length scale, the Kuhn length. Once this length has been rescaled we must be left with a universal model, and 
since the length rescaling in A e /j, n occurs completely independently of the amplitude, this means that A c must be a 
universal quantity. 

We also observe that we now have an exact critical exponent relation, 7 = 1, to which there should be a 
corresponding relation involving amplitudes. 

It would be highly desirable to identify the critical amplitude with a standard quantity in the scaling theory of 
polymers. 
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We note in passing that this universal amplitude may serve as a signal of the 6 transition of interacting self- 
avoiding walks. We know that the amplitude of endless simple random walks is unity. For two-dimensional 9- 
polymers we have 7 = 8/7 [45], and so we could expect that A c for endless ^-polymers would have a distinct value, 
i.e. it should lie between 1 and 1.57- • • . The upper critical dimension for the 9 transition is d = 3, and so for 
three-dimensional endless ^-polymers it is likely that A e = 1. 

We will just make one further point. Amplitudes for lattice models, as opposed to amplitude ratios, are not typi- 
cally universal quantities. For comparison, the amplitudes for SAW are A = 1.1449355(5) [46] for the honeycomb 
lattice, and A = 1.1770423(1) [37] for the square lattice. In contrast, the estimates for A e from the square and 
honeycomb lattices agree to the fourth decimal digit. 

5.2 Corrections to scaling 

The analysis for eSAW on two-dimensional lattices in Sec. 4.5 supports our contention that there are no analytic 
corrections to scaling for the number growth, by confirming that the amplitude ai from (2.1) is indistinguishable 
from zero for all lattices considered. 

With longer enumerations from the finite lattice method it may be possible to estimate a 2 or even a 3 , which could 
strengthen this confirmation significantly. Monte Carlo simulations could also be used to improve the accuracy of 
the estimation of a\. 

We find clear evidence that the leading non-analytic correction to scaling for two-dimensional eSAW has expo- 
nent —3/2, which is the same as for two-dimensional SAW [8]. In particular, we do not see any evidence of the 
correction to scaling with exponent —11/16, predicted by CFT and observed for another recently invented SAW 
variant [12, 13]. 

For three-dimensional eSAW it is certainly possible that either the leading analytic correction to scaling has 
vanishing amplitude, or that there are no analytic corrections to scaling at all. However, it will be far more difficult 
to gather numerical evidence to resolve this question. This is because the leading non-analytic correction to scaling 
dominates the analytic correction, and in addition there are two non-analytic corrections which would compete with 
the leading analytic correction. The leading analytic correction has exponent —1, while the next-to-leading non- 
analytic corrections are — 2Ai ~ —1.06 and — A 2 which is believed to be close to —1. Consequently we consider 
the prospects of obtaining sufficiently strong numerical evidence to determine the amplitude of the leading analytic 
correction to be remote for three-dimensional eSAW. 

Nonetheless, in the absence of a compelling argument for the existence of confluent logarithmic corrections, we 
consider it likely that there are no analytic corrections to scaling for three-dimensional eSAW. 

5.3 Anti-ferromagnetic singularity 

We have examined our series on the bipartite lattices - square, simple cubic, and bcc - for any evidence of the 
anti-ferromagnetic (AF) singularity (the honeycomb lattice is bipartite, but we are only considering one of the sub 
lattices). The AF singularity arises from polygon-like configurations: the odd-even effect is due to the fact that 
self-intersections can only occur between either two odd or two even sites of the walk. 

The signature of the AF singularity for SAW is unambiguous for the same series lengths as we consider here 
for eSAW. From differential approximant analyses one can clearly detect the singularity on the negative real axis, at 
the correct location. From direct fitting, when we include an AF term in our fits we get clear signs of convergence 
to non-zero values for the fitted amplitudes. Finally, the inverses of the generating functions of interest provides a 
strong signal for the existence of the AF singularity [47]. 

For the case of eSAW, none of these analyses are conclusive! 

For the bipartite lattice series, there is no clear signal from the DA analyses. On occasion, the differential 
approximants have singularities on the negative real axis, but the typical behavior is a strong signal of the critical 
point on the positive axis, and then other singularities scattered in the complex plane. 

From the series themselves there are clearly oscillations, but the period in no case appears to be strictly two, and 
we have not been able to usefully apply the direct fitting approach. 

Calculating the inverse of the generating functions, a la [47], leads to a strictly alternating series for SAW, but 
for eSAW on the square lattice there is no regular pattern at all, while for eSAW on the simple cubic and bcc lattices 
there is some hint that the series are settling down to something akin to alternation, but this is by no means definitive. 
Indeed, eSAW on the triangular lattice - which certainly does not have an AF singularity - displays a similar pattern 
to the simple cubic and bcc series. 

In summary, the series analysis is inconclusive with respect to the existence of the AF singularity, and it may 
well be that the observed oscillations are determined by singularities in the complex plane. 
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We will now attempt to explain our intuition for whether there should be an AF singularity based on physical 
arguments. This is purely speculative, and in the absence of further numerical evidence should perhaps be given 
little weight. 

For SAW, polygon-like configurations are the origin of the AF singularity: the dominant contributing configura- 
tions consist of large loops, and the odd-even effect arises from the interaction of the end-points. 

For eSAW, polygon-like configurations are exponentially rare, as explained in Sec. 5.4. Large loops are possible, 
but the ends are then forced apart by the endless condition. Heuristically, for an eSAW of n steps, 0(n) steps can 
be part of a large loop, but then 0(n v ) steps are required for the end-points to escape from each other. We expect 
that this will result in the odd-even effect from the end-point interaction being exponentially suppressed or even 
eliminated completely, and hence no AF singularity will lie on the radius of convergence. It may be that our 
intuition is incorrect and that merely the amplitude of the AF singularity is suppressed, but our best guess is that the 
AF singularity is completely absent for eSAW. 

In future, we will revisit this question with longer series in hand in an attempt to resolve this question definitively. 

5.4 Open questions 

One important topic for eSAW which we have not yet addressed is the behavior of eSAW with fixed end-points, and 
in particular the asymptotic behavior of e n (x) as n approaches infinity with x fixed. For SAW, it is well known that 

c n (x) ~ B(x)n n n 2 - a , (5.1) 

but the situation for eSAW is radically different. For eSAW in d dimensions, in the limit of large n it seems that 
the eSAW must become quasi (d — 1) dimensional objects, with the eSAW effectively confined to a 'slab' with 
dimensions fixed by the end-to-end vector x. This slab could have quite a complicated structure, with the parallel 
surfaces allowed to have arbitrary shape provided that they are always separated by x. Thus it is not clear if the 
resulting eSAW would have growth constant of d — 1 dimensional SAW or a value that is intermediate between 
d — 1 and d dimensional SAW. The challenge is to understand and possibly predict the values of k, a, and B e in 

e n (x) ~ B c {x)n n n a . (5.2) 

If eSAW in this limit are exactly like fixed end-point SAW in d— 1 dimensions, then we would expect that k = fid-i 
and a = —(d— l)z^-i- This is a question that might be resolvable with a good theoretical argument, and it is also 
clearly desirable to study this question numerically. 

Another important question is whether the pivot algorithm [48, 49] is ergodic for eSAW. One reason for the 
importance of this question is that recent improvements to the implementation of the pivot algorithm [50, 44, 51] 
have made it even more powerful. In addition, it is now clear that the pivot algorithm can be used for extremely 
accurate calculations of /j,, 7, and v for SAW. 

We expect that the pivot algorithm is indeed ergodic, since eSAW are very similar in nature to SAW, and in 
some sense they allow for less pathological behavior with respect to the pivot algorithm than SAW. This is true in 
the following sense: for SAW it is possible to form spiral structures with their ends which can only be unfurled 
extremely slowly with pivot moves. In contrast, for eSAW it is never possible for an end to become trapped, and 
thus it would seem that 'most' sites of a typical eS AW should be candidates for valid pivot moves. The key question 
is whether there are novel equivalence classes of eSAW which are invariant under pivot moves, which would result 
in non-ergodicity. 

It may be difficult to prove that the pivot algorithm is indeed ergodic for eSAW. The difficulty lies in the fact 
that the non-local interactions change when a pivot is made, as the end-to-end vector typically changes after a 
successful move. Combining fixed end-point pivot moves [52, 53] with standard pivot moves might allow for a 
more straightforward proof, but this is a sub-optimal solution because fixed end-point moves are less efficient for 
the purpose of simulation. Application of the pivot algorithm to eSAW has not been the focus here, but it is a topic 
we intend to explore in the future. 

Finally, much remains to be done regarding the application of theoretical techniques to eSAW. In particular, the 
renormalization group in conformation space [17] could be adapted to the study of eSAW, and should provide a 
strong theoretical argument that 7 C = 1 in all dimensions. 

5.5 Applications 

Perhaps the application for which eSAW has the clearest advantage over the competing models of SAP and SAW is 
in the study of polymer knotting. Simulations for eSAW promise to be far more rapid than for SAP, provided the 
pivot algorithm with free end-point pivot moves can be proved to be ergodic (see Sec. 5.4). 
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At first glance, the eSAW series seem to hold less information than the corresponding SAW series of the same 
lengths, but it is quite possible that this is a short length effect, and that eSAW becomes superior to SAW at longer 
lengths. The absence of analytic corrections-to-scaling, and the possibility of deriving a theoretical prediction for 
the universal amplitude A e , are significant advantages for eSAW over SAW. The large n limit can be probed via 
finite lattice method enumerations (perhaps approaching n = 80 or so for the square lattice), and also via Monte 
Carlo simulation using the pivot algorithm (where n > 10 8 ). Thus we believe that eSAW hold considerable promise 
for allowing extremely accurate calculations for critical exponents and the connective constant for polymers in two 
and three dimensions. 

Finally, we speculate that eSAW may be useful in the context of exact results on the honeycomb lattice [7, 29]. 
eSAW are a naturally defined special subset of SAW which have zero winding angle. 

5.6 Generalizations 

The concatenation idea of eSAW can be straightforwardly applied to other models of linear polymers, including 
interacting self-avoiding walks (IS AW) [54], the Domb-Joyce model of weakly interacting walks [55], and indeed 
continuum models such as bead-spring models. For these models, the interaction energy would need to be converted 
to an interaction energy per unit length, but other than that there are no significant complications. The self-avoiding 
trail model [56, 57, 58] and generalizations (e.g. [59]) could also be easily made 'endless'. 

One interesting consequence of the existence of a universal amplitude for the endless Domb-Joyce model, is 
that as the interaction parameter is varied, the (weighted) growth constant will vary while the amplitude will remain 
fixed, i.e. we would expect that for interaction strength J: 

e„(J) ~ A el i{J) n . (5.3) 

Of course this exactly what happens for /it( J) and 7 for the usual Domb-Joyce model. 

It would be interesting to see if the eSAW idea could be extended to the continuum Edwards model [60]. 

As flagged in Sec. 2.7, it may be possible to allow for eSAW with an odd number of steps on the honeycomb 
lattice. More generally, the rule for concatenation of a walk with itself could be expanded to allow for symmetry 
operations such as reflections to be applied. 

One weakness of the eSAW model is that there is no natural way to incrementally grow an eSAW. For a SAW, 
any sub-walk is also a SAW, but this is not true of eSAW. Interestingly, the set of eSAW for which any sub-walk 
is also an eSAW is precisely the set of directed walks. This can be seen straightforwardly by noting that: (a) If 
an eSAW has two steps in opposite directions, then the sub-walk which begins and ends with these steps is not an 
eSAW, and (b) Every directed walk is an eSAW. Perhaps there are other subsets of eSAW which may be of more 
interest. 

It may be possible to study polymer networks using variations of eSAW, by concatenating together repeated 
copies of star or branched polymers. An example of a two-dimensional network created by such a recipe is shown 
in Fig. 5.1. 




Figure 5.1: An example of a self-avoiding network created via concatenation of a star polymer on the square lattice. 
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6 Conclusion 



We have introduced a new model of self-avoiding walks for which we have theoretical and numerical evidence that 
end-effects should be completely absent. 

The enumeration results and analysis confirmed our expectations that 7 = 1 for eSAW, but more serious 
numerical work should be done to provide stronger confirmation. In future we will apply Monte Carlo tech- 
niques [49, 50, 44, 51] to accurately estimate critical exponents. It would be desirable to improve the rigor of 
our heuristic arguments that 7 e = 1 and u c = v. 

We have made the surprising discovery that the eSAW model has a universal amplitude, and estimate that A c = 
1.57075(10) in two dimensions, and A c = 1.183(3) in three dimensions. We have given one argument for this 
phenomenon, but it would be highly desirable to find other theoretical justifications. It remains to be seen if other 
models can be found which possess universal amplitudes. We also believe that there are no analytic corrections to 
scaling for the number growth of eSAW, and have supporting numerical evidence for two-dimensional eSAW. 

The concatenation technique described in this paper can be straightforwardly applied to other models of poly- 
mers to create new 'endless' models. We are hopeful that the simple idea introduced here will open up a fruitful 
avenue of research. 

A Enumeration data 
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16 
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18 
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114635882508 


1 106794584 


26352252 


22 
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53595401400 
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27 
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6906918040 


28 
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2595592842456 


44751600732 


30 


6847276871420 


1213539966457200 


6847273559040 


114121225984 


31 


18063481834652 


3363306014128316 


18063481834648 


291346481204 


32 


47652479932748 


9306731807696384 


47652471265792 


744569863528 


33 


125710536151052 


25715129922994380 


125710536083088 


1904705092168 


34 


331633319500804 


70954363887188496 


331633296639192 


4876960244694 


35 


874875502020444 


195525476447120540 


874875502018840 


12498221457412 



Table 2: Enumeration data for the square lattice. 
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74100 


10 


6235950 
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3312601632 
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17 
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18 


1429710254742 
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1429708907808 
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19 


6693530115990 
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176145529368 


20 


31322452467006 


1417625514526560 


31322446230480 


783061155762 


21 


146644481901582 


7039475232292062 


146644481840472 


3491535281916 


22 


686366924272686 


34860099689171064 


686366895112272 


15599247616188 


23 


3213444320539710 


172198320473193822 


3213444320539704 


69857485229124 


24 


15042778398965358 


848672945442198144 


15042778262399904 


313391213799998 



Table 3: Enumeration data for the simple cubic lattice. 
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12 


2 


18 


30 


56 


56 


132 


3 


60 


114 


296 


344 


1284 


4 
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486 
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12756 


5 


726 


1986 
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14648 


126612 


6 


2448 


8202 


62648 
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1265604 


7 
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33858 


365464 
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8 
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140262 


2130488 


4066376 
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9 
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580938 
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26456024 


1267690812 


10 


336618 
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72464136 


172470776 
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11 


1150320 


9980394 
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1124590904 


127218134748 


12 
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41384106 
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7328945336 


1274923451364 


13 


13419204 


171634482 


14384367304 


47829176312 


12779117065740 


14 


45828192 


711932538 
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311758824440 
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15 


156512220 


2953394214 


489611186736 


2035191240584 




16 


534463698 


12253006566 


2856635570680 


13269735071816 




17 


1825120584 


50838740046 


16667509269248 


86634355605416 




18 


6232259412 


210946278126 


97251708586688 


565037351668904 




19 


21281168202 


875324356398 








20 


72666555570 


3632308990686 








21 


248124503652 


15073375131306 








22 


847224827676 


62553271396074 








23 


2892836367066 










24 


9877456541376 










25 


33725891989626 











Table 4: Enumeration data for the honeycomb, triangular, union jack, bcc, and fee lattices. 
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